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f^ 1 We investigate the Lorentz contraction of a deuteron in flight. Our starting point is the 

Blankenbecler-Sugar projection of the Bethe-Salpeter equation to a 3-dimensional quasi 
potential equation and apply it for the deuteron bound in an harmonic oscillator potential 
(for an analytical result) or by the Bonn NN potential for a more realistic estimate. 
We find substantial quenching with increasing external momentum and a significant 
modification of the high momentum spectrum of the deuteron. 
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1. Introduction 

One of the main goals of physics with modern accelerators is the detailed investiga- 
tion of the response of the nuclei^ to large momentum transfers of around or above 
1 GeV/c . In particular, here the lightest nucleus, i.e. the deuteron, is a preferred 
object due to its well known nuclear structure . In recent years many investigations 
with the deuteron as a target have been performed at large momentum transfers 
with proton and electron beams^, further studies with deuteron beams have beerPl 
or will be performed^ in the near future. Here the deuteron as a spin 1 object is 
very special, as it allows to measure in addition to the total and differential cross 
section and the analyzing power also the tensor polarization of deuteron beamg^. 
It is hoped - and that explains the genuine interest - that high momentum transfer 
processes reveal details of the short range reaction mechanism and of the deuteron 
itself, which are expected at momentum transfers up to the low GeV/c region to be 
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dominated by meson exchange currents" and the excitation of baryon resonancea^. 

In most calculations the deuteron is conventionally described as an object in rest 
with standard s and d wave components. However it is clear that with increasing 
momenta or momentum transfers on the deuteron this picture has to be modified: 
in coordinate space a moving deuteron is quenched along its axis of motion. As a 
consequence the deformation of the deuteron is increasing with increasing momenta 
which leads to a significant modification and enhancement of the deuteron spec- 
trum at high momentum components. This, as expected, will modify substantially 
predictions of observables, which depend sensitively on large momenta components 
of the deuteron. 

In standard calculations the quenching of the deuteron is frequently not in- 
cluded, instead one uses approximate descriptions, such as to go to the Breit frame 
for the deuteron form factoJJJJ, in order to hopefully minimize the error due to the 
approximation. 

In this note we investigate the quenching of a moving deuteron in a simple model. 
In chapter 2 we start from the manifest covariant Bethe-Salpeter equationAA; and 
solve, after projecting onto 3 dimensions, a still covariant quasi potential equation^ 
for an harmonic oscillator potential of the deuteron with the goal to obtain analytical 
solutions. These solutions are then derived explicitly within a few approximations for 
the d-state of the deuteron, and presented in chapter 3, together with characteristic 
results for the deuteron wave function in momentum space; subsequently we apply 
the same prescription to a realistic deuteron wave function from the Bonn potential. 
Finally we close in chapter 4 with a summary and some remarks to improve of the 
approach. 

2. Basic equations 

In this chapter we derive for a given kinematics the quasi potential equation for 
the deuteron from a Blankenbecler-Sugar projection of the Bethe-Salpeter equation 
(BSE). 

2.1. Bethe-Salpeter equation for the deuteron 

For the deuteron with the total momentum P and the relative momentum q between 
the two nucleons the BSE is given as 

(£ + ^ + M)(£_^ + M) f 

thiP, q) = —. ^^ ^ -^T^ / K{P, q, kU(P, k)dk (1) 

where K{P, q, k) denotes the kernel of the BSE and M the nucleon mass. The total 3- 
momentum P is specified along the z-axis, i.e. P = (0, 0, Pz), while the 4-momcntum 
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P is given hy P — {Ed,0,0,Pz) with the total energy of the deuteron given as 
Ed = \J M'd + -P| (Md is the deuteron mass). To simphfy the comphcated spinor 
structure of '0(^; <?); which is a 16 component spinor, we neglect in M^ = 2M + e^ 
and 2M ^ Ed the binding energy of the deuteron Ed = —2.2 MeV, and put the two 
nucleons on the mass shell; then (-j ± ^ + M) = 2M (this corresponds in keeping 
only the large components in each of the Dirac spinors of the two nucleons). With 
this approximation the BSE is simplified drastically to 

UP 



ii^ + q)2-AP-ivm-qr-M^-ir,) 

X / K{P, q, k)ij{P, k)dk . (2) 



2.2. The Blankenhecler- Sugar (BS) projection 

As the BSE can - except for separable potentials and a few academic cases - only 
be solved numerically, we simplify the problem by a reduction of the BSE to a 3- 
dimensional, but still covariant quasi-potential equation (QPE) by restricting the 
nonstatic effects in the BSE. Unfortunately the transition from 4 to 3 dimensions 
is not unique: there exists an infinite number of QPIji^. Thus we have to invoke 
- guided by the underlying physics ~ plausible arguments for the reduction of the 
BSE. As the two nucleons in the deuteron are expected to exhibit a symmetrical de- 
pendence on the fourth component of internal momentum g, we put them, following 
the arguments of Blankenhecler and Sugail^, in their Greens function 

AIVP 
^^^' '^^ " ((£ + g)2 - M2 - i,y) ((f - g)2 - M2 - ir^) ^^^ 

symmetrically off their mass shell 



((f + qY - M2 - ir;) + ((f - g)2 - M^ - ir?) 
x^Q((f +9)^-M^)-((f-9)^-Af^)) (4) 

where the (5-function can be solved for the go-component to give 

5{Pq) ^ ^5{q, ^q.). (5) 

As a consequence the Greens function is now reduced to 

Ggp{P, q) = 2i7r— — -5 i^ ^ (6) 



Er 



{^-mr<i'^-i. 



KP 



which depends explicitly on the external momentum Pz of the deuteron. We verify 
that in the limit of a deuteron in rest, i.e. Ed — > Md — 2M -f Ed the denominator D 
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in the Greens function recovers the standard non relativistic structure 



\ 4 



In this picture the full deuteron Greens function in the QP limit is then given 
dropping terms quadratic in Sd, as 

G{P,q) = ■nr——6{qo - -^qz)— m,^2 2 JT 

21^d i^d [Med - (-57)9? - 1i ) 



. Md SJlo-^lz 

Ed\Med-{%^Yql~q\y 



mMC-^) ''\, ^''''" . (8) 



With these steps the structure of the QPE is, for a given interaction kernel 
K{P,q,k), completely specified (the final expression for the QPE is given in the 
following chapter 3). 

3. Momentum spectrum of a deuteron in flight 

In this chapter we derive the structure of the momentum dependent harmonic kernel 
and present the solution of the corresponding QPE both for the central and the 
tensor piece of the deuteron wave function in momentum space. 

3.1. Harmonic kernel in momentum space 

As already mentioned we use for the potential of the deuteron an harmonic kernel, 
in order to end up with analytical solutions of the QPE. Thus we start from the 
four-dimensional form 

K{x) = V{x) - Xho^'' (9) 

or its Fourier transform 

Viq -k)^ Xho /"a;2e'(«-'^-)^da; 

= -i2n)^XHo-^Siq~k). (10) 

To reduce the 4 dimensional ^-function to a 3-dimensional form we start from the 
representatiorii2. 

dq^ m^o d'm? (q — k^ — "m? — irj 

or, upon performing the BS-projection 

-—S{q-k) = -l[ia 4-^^. (12) 

d«2 ^ ' m-.0am3(ML)2(^ fc )2 + ( ;, )2+^2 ^ ' 
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so that the harmonic interaction kernel V{q — k) reads with the notation 

Md Md 



Ea yWI+Pf) 



A (13) 



V{q-k) = {2nfXHO^(^^ + ^]s{q-k). (14) 



Md \ A2 rfg2 aq2 



3.2. Deuteron wave function 



For a deuteron in rest the total wave function of the deuteron is given schematically 
as a combination of s and d wave pieces 

-ipiEd^q) = Csipsiq) + CdSuigjipdiq))- (15) 

Thereby c^ and Cd define the strength of the s and d-wave, which depend on the 
underlying nucleon-nucleon potential; with the tensor operator 
'S'i2('i') == 3o'i(70'29 ^ o'iO'2, this yields in a coupled form 

^{Ed,q) = (c,i?,(g)[yo(i)[H]']'" + CdRd{q)[Y2{qW^VV"') [ii]°'° (16) 

with radial functions Rs^d{q) (the coupled form [^i]'^''^ denotes the isospin content 
of the deuteron) . In the following we discuss the solution for the s and the d-wave 
separately. 



3.3. Solution for the s-state 

With the steps given above, we obtain from Eqs. (8,15) and with the additional 
notation 

A2 = i(2^)5AHo (17) 

the final equation for the s-wavc in momentum space 



dq\^ Ed 



TTTT + A't-T - Msd - {-^fql - ql ^P, q) = 0. (18) 



Note that the notion s-wave is somewhat misleading, as the motion of the deuteron 
along the z-axis breaks the rotational symmetry of the standard s-wave component 
and introduces angular momenta I ^ 0. 

The solution of Eq. (18) is easily obtained, following standard text books. Going 
over to cartesian coordinates with 

:^^+A^^+A^^)-Me,A2g|-g^-g^))^(P,g)=0 (19) 
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we obtain with the notation 



^z- J2^ ^rr - ^ - \ (20) 

Sd = ex+Sy+ Ez (22) 



the final equations 



^ + fc.'-Ak')0(9.)=O (23) 

^ + kf-XJqf^ct>{q,)=.0 (24) 

(where i — x,y). Then the solution for the radial part is given as 

(/.(ft) = 7V,e5^'«- (25) 

(Ni are normalization constants). 

The modification of the momentum spectrum along the z-axis is easily seen in 
comparison with a deuteron in rest: comparing with the solution of the spherical 
harmonic oscillator 

</. ~ e"^ (26) 

where a is the common oscillator parameter, we find from 

K--\y = j = a^ (27) 



^^(Md,2 2 

A ^Ed 



A, - — - (i^)'a' < \cc,y forP, > (28) 



immediately an enhancement of the z-components of a deuteron in flight or, corre- 
spondingly, a quenching of the z-component of the deuteron in coordinate space. 

3.4. Solution for the d-state 

The d-state component of a deuteron in flight is more complex than the correspond- 
ing s-statc. A rigorous solution, even for an harmonic kernel, can be obtained only 
numerically, as the tensor operator couples longitudinal and perpendicular momen- 
tum components. We suppress this coupling by the assumption that in 

Si2{q} = V2l^{ql + ql)[Y2{qz4j[a„a,]Y° (29) 

for \ q \^ qz ^ the case we are most interested in - the unit vector for the total 

momentum is directed along the z-axis (with Y20 — \/ ^ a-nd F22 — v i^) then 
only the terms q"^ and q^ without an explicit dependence on the angular momentum 
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enter in the decoupled radial equations for q^ and q . Consequently we end up, after 
separating into cartesian coordinates, with the radial equations 

iq■-f2+kf->^hf)Hq^) = o (30) 

dqf 

for the x,y and z components (following the notation for the s-state equation in 
3.3). The solutions of these differential equations are readily obtained: the general 
equation 

x\f"{x) + (fe2 - a^x^)f{x) = (31) 

is reduced with suitable substitutions 



f{x) — ^/xy{x), u = y —d?-x (32) 

to the differential equation for Bessel functions 

u^y"{u) + uy'{u) + [v? -{\- \?))v{u) = (33) 

with the corresponding regular solution 

/(x) = V^J(i_,,^ (y^x) . (34) 

(J is a Bessel function of first kind.) 

3.5. Results and discussion 

In this section we present typical results from the formalism derived in the proceed- 
ing sections. Thereby we focus on the momentum spectrum of the s and d-state, as 
here the influence of the quenching of a deuteron in flight is seen most clearly. To 
model our static deuteron we use an oscillator parameter a — 0.7 fm as a typical 
scale (of course we are aware that the representation of deuteron wave function as an 
harmonic oscillator with a single size parameter is an inadequate representation for 
a realistic comparison). Typical results are given in Fig. l(a,b) for the s and d-wave 
momentum spectrum respectively. We find that at momenta around IGeV/c for a 
deuteron with a longitudinal momentum of IGeV/c - a typical momentum antici- 
pated for the deuteron beam of COSY ~, the spectrum both in the s and d-wave is 
enhanced by typically one order of magnitude. In addition, the spherical symmetry 
in the s-state is increasingly broken with increasing P^ . This is most readily seen in 
coordinate space, where one obtains upon the corresponding Fourier transform 

mJ + pJ 



M^,r)^e'-'\ -T-'+'-y (35) 



or, if we introduce explicitly angular momentum functions, 



(36) 

Clearly, as easily seen from an expansion of the exponential function, for Pj > the 
s-state wave function for a deuteron contains in principle an infinite sum of even 
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Momentum Spectrum of Deuteron, Harmonic Potential, S - State 



P = OGeV/c ■ 
1 GeV/c 
2GeV/c ■ 




12 3 4 5 
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Fig. la. Momentum spectrum of the s-wave of the deuteron for different external momenta Pz 



P = GeV/c ■ 

1 GeV/c 

2 GeV/c ■ 




Fig. lb. As Fig. la, however for the d-wave of the deuteron. 



angular momentum components; spin-orbital coupling admixes an additional L=2 
component for the s-state wave function. A similar qualitative argument holds also 
for the d-statc (quantitatively, it differs by the additional coupling of 12™ (^) from 
the structure of the d-state in the static deuteron). We extent our results for an 
harmonic kernel to deuteron wave functions from modern NN potentials, using the 
recipe q^ — > X'^{Pz)q^ explicitly for the deuteron wave function from the Bonn po- 
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Momentum Spectrum of Deuteron, Bonn Potential, S - State 




Fig. 2a. Momentum spectrum of the s-wave of the deuteron for different external momenta Pz 
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Fig. 2b. As Fig. 2a, however for the d-wavc of the deuteron. 



tential. The result is shown in Fig. 2(a,b) for the s and the d-state of the deuteron. 
QuaHtatively, the resulting pattern - the enhancement of the Qz spectrum for P^ > 
- is similar as for an harmonic kernel, though quantitatively such a modification is 
less pronounced, as already for a deuteron in rest the Bonn potential significantly 
enhances the high momentum spectrum of the deuteron wave function relative to 
an harmonic potential. 
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4. Summary, conclusions and outlook 

In this note we investigate the momentum spectrum of a deuteron in flight. Starting 
from the BSE and foUowing the procedure of Blankenbecler and Sugar, we derive a 
3-dimensional QPE for the deuteron. Introducing an harmonic kernel for the bind- 
ing of the deuteron, we solve the corresponding equation for the s and d-state of 
the deuteron analytically. 

As the main result we find, that for the longitudinal component (along the exter- 
nal momentum Pz) the momentum spectrum, if compared with a deuteron in rest 
or with the perpendicular components, is substantially modified with increasing P^: 
the high momentum components are significantly enhanced, as expected from the 
quenching of the deuteron in coordinate space along the direction of motion. This 
is not surprising, as even in the s-state the spherical symmetry is broken for P > 
and higher partial waves I > modify the s-state (with / = for a deuteron in rest), 
whereas the conventional d-state develops in addition an 1=0 component. Thus the 
deuteron gets more and more deformed with increasing P^. 

From the quantitative point of view our results are certainly preliminary: clearly 
our approach has to be based on a more realistic interaction kernel for comparison 
with a deuteron in rest, by taking the deuteron wave function or deuteron form 
factors from modern nucleon-nucleon potentials^! A simple, but still approximate 
description is to fit for example the elastic deuteron form factor for a deuteron 
in rest by a set of Gaussians with different size parameters along the results in 
chapter 3 (the result is then still analytical) . A more rigorous step would be to solve 
d-quenching for a kernel from modern nucleon-nucleon potentials, at the expense 
however, to end up with a complicated system of coupled differential equations due 
to the coupling of longitudinal and perpendicular components for Pz > 0, which 
can be solved only numerically, so that the transparency of the analytical solution 
is lost. In addition it has to be stressed that Lorentz boosts in the instant form are 
dynamical and involve the underlying interaction- -^ (opposite to kinematical boots 
on the light cone ) ; thus there exist different prescriptions for their incorporation 
| J-g | J-V I Experimentally, detailed information on Lorentz boosts of the deuteron will 
result a systematic study of pd and, in particular, dp reactions at COSY, where a 
(polarized) deuteron beam will be available in the very near future 
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